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Abstract 

In  order  to  get  the  corresponding  relationship  of  target  distance  and  the  initial  speed  of  muzzle  at  the  moment  of  non- 
lethal  blow  of  constant  terminal  effect,  a  certain  type  of  9mm  pillar  rubber  kinetic  energy  ammunition  is  taken  as  an 
example.  The  Mutilating  physiological  effect  experiment  and  numerical  fitting  are  conducted  to  comprehensively 
analyze  the  Efficiency  value  when  being  struck.  With  exterior  trajectory  equation,  compute  the  exterior  trajectory 
according  to  Siacci  Method  and  validate  with  simulation.  Then  polynomial  fit  modeling  is  applied  to  the  calculated 
statistics  to  get  the  simple  empirical  equation  corresponding  to  the  target  distance  and  the  initial  speed  of  the  muzzle. 
Based  on  this,  the  father  interior  ballistics  design  of  non-lethal  blow  of  constant  terminal  effect  and  energy  control 
experiment  of  kinetic  energy  ammunition  could  be  well  done,  and  the  exterior  trajectory  models  and  empirical 
equations  of  other  types  of  kinetic  energy  ammunitions  could  be  calculated. 
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1.  Introduction 

When  police  executes  anti-terrorism  task,  main  rioters  appear  inconstantly  and  the  location  is  changing 
continuously.  In  order  to  effectively  combat  with  rioters,  control  development  of  the  situation,  special 
requirements  on  weapon  equipment  are  proposed.  Killing  power  is  adjustable,  and  terminal  performance  is 
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controllable.  The  current  improved  design  is  on  basis  of  existing  domestic  riot  gun.  According  to  the 
distance  of  target,  opening  in  the  appropriate  parts  of  barrel  is  to  vent  out  right  amount  of  gunpowder  gas, 
adjusting  kinetic  energy  projectile  muzzle  velocity  is  to  achieve  constant  terminal  efficiency  non  fatal  blow 
of  different  distant  target. 

The  basis  of  adjusting  kinetic  energy  projectile  muzzle  velocity  is  accurately  obtained  the  relationship 
between  target  distance  and  muzzle  velocity.  The  paper  takes  a  9mm  cylindrical  rubber  kinetic  energy 
projectile  for  example,  and  establishes  outside  trajectory  calculation  mathematical  model  to  deduce 
empirical  formula  about  the  relationship  between  target  distance  and  muzzle  velocity,  which  can  lay  the 
foundation  for  internal  ballistic  analysis. 

1.  Determine  the  constant  terminal  performance  value 

Riot  kinetic  energy  weapon  combat  against  staff  who  expose  skin  and  wear  a  small  amount  of  clothes, 
projectile  directly  effect  on  target  skin,  universal  threshold  can  be  extracted  with  kinetic  energy  for  the 
standard.  It  should  adopt  the  specific  kinetic  energy  injury  criterion  [1],  While  striking  target,  the  ratio 
between  kinetic  energy  of  the  projectile  on  the  contact  with  the  maximum  area  of  the  target  sectional 
surface  parallel  to  the  projectile  is  specific  kinetic  energy  (Ed). 

Firstly,  the  damage  physiological  effect  experiment  carries  out.  The  test  uses  biological  targets, 
including  pigs,  leg  of  pork,  pork  chops,  etc.  The  experimental  device  is  the  United  States  miro3  type  high 
speed  camera  and  XGK-2002  small  target  speed  measuring  laser  screen  target.  The  ballistic  gun  is  fixed, 
and  kinetic  energy  projectile  launch  platform  is  established.  By  adjusting  loading  dose  and  different  target 
distance,  the  different  target  speeds  are  received  to  achieve  combating  against  the  layer  skin  of  biological 
targets  in  different  specific  kinetic  energy  [2],  The  analysis  results  are  as  follow  Table  1. 

Table  1.  penetrate  depth  of  the  projectile  of  different  specific  kinetic  energy 


specific  kinetic  energy  <7  8 

Ed  (J/cm2) 

!  9  10 

11  12  13 

penetrate  depth  0  C 

1  0.2  0.6 

1.2  1.9  2.9 

Combined  damage  physiological  effect  experiment  with  foreign  relevant  experimental  results  [3-4],  the 
analysis  is  that  when  the  specific  kinetic  energy  reached  to  10J/cm2  left  or  right,  the  projectile  was  bounced, 
but  began  to  invade  the  skin  of  biological  targets  and  just  broke  it.  The  hit  parts  are  tumid  and  slightly 
damaged,  which  is  consistent  with  mild  damage  standard.  Considering  nonlethal  effects  and  lethality  of  the 
kinetic  energy  projectile,  when  blow  constant  terminal  kinetic  energy  is  10J/cm2,  the  projectile  can  make 
the  target  loss  resistance  instantly  without  causing  greater  harm.  Therefore,  when  the  constant  terminal 
performance  is  10J/cm2,  the  projectile  velocity  is  calculated,  v  =79.740m/s. 

2.  Computation  of  Exterior  Trajectory 

Given  the  range  and  kinetic  energy  of  projectile  are  not  too  high,  the  whole  movement  is  the  movement 
near  the  surface,  do  assumptions:  the  projectile  is  axisymmetric,  the  nutation  angle  8  =  0;  as  the  small  range, 
do  not  take  the  curvature  of  the  earth  and  changes  of  the  acceleration  due  to  gravity  with  height  (g=9.8m/s2, 
and  the  direction  is  always  vertical  down);Coriolis  acceleration  ac  =  0;standard  meteorological  conditions: 
no  wind,  rain  and  snow.  Under  these  assumptions,  the  forces  acting  on  the  projectile  only  are  gravity  and 
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air  resistance,  regard  the  projectile  as  a  particle  whose  mass  being  concentrated  in  the  center.  Get  the  vector 
dv 

equation  of  the  motion:  —  =  a  r+  g  ( 1 ) 

dt 

Project  the  vector  equation  (1)  on  the  horizontal  direction,  then  get  the  exterior  trajectory  equation: 


~r  -  —CH  (y)G(v) 
dx 

(2) 

d0  __  g 

(3) 

dx  v2 

*=!»« 

(4) 

dx 

dt  1 

(5) 

dx  u 

Thereinto, 

(x,y) - coordinate  of  the  projectile  at  any  trajectory  point;  v - the  speed  of  the  projectile;  C - 

Ballistic  coefficient;  H(y) - a  function  of  air  density,  low-stretch  ballistic  standard  conditions:  H(y)=  1; 

G(v) - air  resistance  function;  G(y)=4.737x  \OavCx0^(M)  ,  CxqN(M)  as  the  standard  missile  drag 

coefficient;  g - acceleration  due  to  gravity,  g  =  9.8m/s2;  0 - angle  between  projectile  axis  with  the 

horizontal  at  any  trajectory  point. 

By  Siacci  Method,  any  point  on  the  trajectory  there  is  0<y>  9 ,  the  introduction  of  U(u  = _ - _ ) 

cosd0 

instead  of  v(v  =  u  ), separate  variables  u  with  0  in  G(  u  &  ),  then  U>  v.  For  the  trajectory  with  small 

angle  (  0  0  <5  °  ),  regard  as  cos  7  0~cos  9 ,  then  error  between  G(U)  and  G(v)  is  smaller,  G(U)  is  an 
increasing  function  of  U,  G(v)<G(U),  for  the  introduction  of  close  substitutes,  with  cosOo  to  further  correct 
the  error.  That  G(v)  =cos00  •  G(U),  by  the  equation  (2),  ^  =  ^4.737*1  IT4  CCx0N(M)dx  ,  That 

k=  4.737  x  10-4  CCx0N(M) ,  then  ^  =  -kdx ,  Which,  G0=v0  ,  U  =  v  ,  We  get  v  =  v{)e^  (6), 

U  COSfc'o  COS0 

take  (6)  into  (3),  de  = - —  dx  ,We  get  tan6i  =  ~g^2fa  ~l)  +  tan0  (7),  compare  (6)  and  (7), 

cos2  6  cos 20ovo2  2kv02coi00 

then 

COS^q  _fa 

y  cose  v°e 

[tan  0  =  ^^- - — +  tan  0O 

2 kv02  cos2  0O 

Table  2.  The  desired  muzzle  velocity  corresponding  to  target  distance 


i  target  distance  (  m) 

muzzle  velocity  Vo,-  (  m/s  ) 

1  10  88 

2  20 

97 

3  30 

107 
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4  40  118 

_5 _ 50 _ 130_ 

_6 _ 60 _ 144^ 

J_ _ 70 _ 160_ 

J _ 80 _ \11_ 

9_90_ 195 


For  constant  kinetic  energy  weapons,  when  firing  with  anti-riot  nonfatal  blow,  effective  range  is  no  more 
than  100m,  for  the  9mm  rubber  cylindrical  projectiles,  m= 2g,  va=230m/s,  rf=9mm,  i= 2,  0(j=%l\  70.  v,  d0,  x,  C 
is  constant,  we  can  find  v0.  Calculated  data  are  shown  in  Table  2. 


3.  Trajectory  Simulation 


With  the  exterior  trajectory  equation,  use  integration  modules  when  establishing  differential  equations, 
determine  the  initial  value  according  to  the  initial  conditions.  Establish  simulation  model  with  Math 
Operation  Module  group  and  Integrator  module  of  Continuous  Module  group.Display  the  simulated 
ballistic  curve  and  speed  curve  with  the  distance  with  the  X-Y  Graph  block.  Simulation  model  is 
diagrammed  in  Fig.  1.  (a).  Write  M  function  file,  set  the  simulation  module  parameter  values  with  the 
function  set_param();  execute  Simulink  model  by  sim  command.  Set  the  initial  velocity  from  velocity  vOi 
in  Table  2;  use  adaptive  variable  step  integration  method  of  odel5s  for  simulation;  with  allowed  relative 
error  range  of  le-10;  start  from  time  0;  end  at  the  corresponding  target  distance.  When  velocity  v0;  less  than 
150m/s,  CxON(M)=0.255-,  velocity  v0i  higher  than  150m/s,  Cx0N(M)= 0.259.  Respective  Simulation  curves  of 
speed  and  the  corresponding  range  are  shown  in  Fig.  1.  (b),  Fig.  2.  (a). 

The  projectile  velocity  as  in  Table  2  obtained  v0l,  in  the  corresponding  target  distance,  the  target  speed 
can  achieve  80m/s  with  a  little  smaller  error.  The  simulation  results  can  reach  the  desired  requirements,  and 
comparing  with  theoretical  computation  values,  it  is  similar  to  a  high  degree  of  numerical  results.  Therefore, 
this  external  ballistics  solution  method  and  computational  results  are  correct. 


llz 

xs 

+ 

1+4 

■■ 

4.  Polynomial  Fitting  Model 

For  the  data  from  computational  results  (xt,  v0i)  (7=1 ,2,. ..,9),  find  a  Polynomial  m  (m  <9),  /?nl(x)=  ao  T 
aix  +  ...  amxm,  so  2>/=2>.-o  m(Xi)Y  =  f(a0,-" am )  (8)  minimum  that  is,  select  at  ( i  =  0,1 

make  f(a0,a/,---am)  =  min  'V  [v,  -  (Hx,)]2  (9)>  where  H  is  the  group  of  polynomials  up  to  m  times.  Satisfy 

y,sH  x-j 
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v. - ^ akxtk jx/  =  0 (/=0, 1 , ‘ •  m)  transpose,  ^ ak xf +y )  =  ^ vtx{  (/=0,l,"-w)  (10).  As  a 

result,  data  processing,  model  fitting,  polynomial  curve  fitting  obtained  as  Fig.  2.(b), 
j7i6)=1.333jd-68.444  (11) 

j726)=0.007jc2+0.628;d-81.381  (12) 

It  can  be  seen  from  curve  of  fitting,  the  result  of  polynomial  fitting  3  or  4  times  is  coincident  with  fitting 
2  times.  Considering  the  fitting  accuracy  and  computing  complexity,  the  kinetic  energy  for  the  9mm 
projectiles,  external  ballistics  modeling  take  m= 2  for  the  best.  For  j72(jc), verify  the  calculated: 
e.«<=|P2(xi)-v0i|  <  Q.5%,  (M,2,...9) 

Accuracy  of  data  fitting  model  meet  the  actual  demand,  so  the  empirical  formula  about  distance  and 
muzzle  velocity:  p2(x)=0.007x2+0.628x+8 1.381 


Fig.  2.  (a)  simulation  curves  (vft>150m/s);(b)  curve  of  polynomial  fitting 


5.  Conclusions 

(1)  Get  the  corresponding  relationship  of  target  distance  and  the  initial  speed  of  muzzle  at  the  moment  of 
non-lethal  blow  of  constant  terminal  effect.  Based  on  this,  the  father  interior  ballistics  design  of  non-lethal 
blow  of  constant  terminal  effect  and  energy  control  experiment  of  kinetic  energy  ammunition  could  be  well 
done. 

(2)  This  method  for  modeling  on  other  non-lethal  anti-riot  low  kinetic  energy  projectiles  is  generally 
applicable.  There  is  a  larger  value  for  the  computation  of  exterior  trajectory  models  of  other  types  of  kinetic 
energy  ammunitions. 
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